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Abstract. In this paper we consider the existence problem of cubature formulas of degree 4fc + 1 
for spherically symmetric integrals for which the equality holds in the Moller lower bound. We prove 
that for sufficiently large dimensional minimal formulas, any two distinct points on some concentric 
sphere have inner products all of which are rational numbers. By applying this result we prove that 
for any d > 3 there exist no a!-dimcnsional minimal formulas of degrees 13 and 21 for some special 
integral. 
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1. Introduction. A main problem of numerical integration is to approximate 
the integral 



Here a; is a d-dimensional vector and /i is a positive measure on a domain f2 in K d . 
We search for an approximant by taking a positive linear combination of the function 
values of / as specified points X = {x±, . . . , xn}, that is, 



We call a cubature formula. The values Wi are the weights and Xi are the points 
of a cubature formula. From the practical viewpoint of numerical computation, it is 
desirable that every coefficient is positive; see [26 , p. 12]. To each formula we assign the 
set of functions for which it is exact. Most often this set is the space of all polynomials 
of degree at most t\ in this case a cubature formula is said to be of degree t. We refer 
the readers to the comprehensive monographs |T0l [24] for the basic theory of cubature 
formula. Hereafter we only consider cubature formulas of odd degree. 

It is well known [53] that the smallest possible number of points |X| in a cubature 
formula of degree 2k + 1 is bounded from below: 
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\X\ > dimP fe (fi). 



(1.2) 
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Here 7^(0) is the space of polynomials of degree at most k restricted to ft. The bound 
(|1.2[) can be improved for some broad class of integrals called spherically symmetric 
integrals [TTIHH] (see Section [5] for the definition of spherical symmetry): 



Here V^ity is the subspace of Pfe(f2) consisting of even polynomials or odd polyno- 
mials if k is even or odd, respectively. A cubature formula is said to be minimal if 
equality holds in one of the above lower bounds. 

A fundamental problem in combinatorics and numerical analysis is the existence 
and nonexistence of minimal formulas for spherically symmetric integrals. In the 
one-dimensional case there always exists a minimal formula of odd degree, and the 
points and the weights of the formula can be uniquely determined in terms of roots of 
an orthogonal polynomial [10 . This one-dimensional formula is known as Gaussian 
quadrature formula. The problem becomes much harder in multidimensional cases. In 
the two-dimensional case many minimal formulas have been already found [U [11] [26l 
|2"TI [3D] ■ An outstanding result is due to Xu [3UJ : He found some special spherically 
(circularly) symmetric integral for which there exists a minimal formula of general 
degree. In higher dimensional cases, however, only a few minimal formulas have been 
found so far, all of which have small degrees or low dimensions. It seems to be the 
conventional belief (see, e.g., E H21 EUJ [21] [32] ) that there exists no minimal formula 
of degree Ak + 1 for any d-dimensional spherically symmetric integrals for k > 1 and 
d > 3 with some possible exceptional cases. There are some reasons why we focus on 
the degree Ak + 1 case: The structure of minimal formulas is strongly restricted. For 
instance a famous theorem by Moller [16] states that the points of a minimal formula 
have the antipodality and contain the origin. 

There are some recent papers involving small degree cases which support the 
conventional belief mentioned above. Hirao and Sawa [T2] observed the relationship 
between minimal formulas of degree 5 and certain combinatorial objects called tight 
spherical 5-designs, and thereby showed nonexistence of minimal formulas of degree 5 
with possible exceptions d = (2m + l) 2 — 2, m > 5; see also [211 [28]. I n the same paper 
they also proved there exists no minimal formula of degree 9 for classical integrals like, 
Gaussian integral on R d , Jacobi integral on the ball. Bannai and Bannai [4 improved 
this result for all spherically symmetric integrals and concluded the conventional belief 
is true for the degree 9 case. However as far as the authors know, very little is known 
on the existence and nonexistence of minimal formulas of degree at least 13 for general 
values of d > 3. 

The present paper has several aims. The primary one is to develop a new theory 
to number-theoretically characterize the structure of minimal formulas of degree 4fc+l 
for spherically symmetric integrals. To do this, we unify a theorem by Mysovskikh [19] 
which is classical in numerical analysis and the Larman-Rogers-Seidel theorem |14) 
(and its generalization by Nozaki |22j). The Larman-Rogers-Seidel theorem is widely 
accepted in combinatorics, but is not fully recognized in numerical analysis and re- 
lated areas. Conversely, the Mysovskikh theory is not fully understood by researchers 
in combinatorics. Thus the secondary aim of this paper is to inform the Larman- 
Rogers-Seidel theorem to researchers in numerical analysis, and conversely, to let 
combinatorialists know details of the Mysovskikh theory. The final aim is to discuss 
the usefulness of our new theory. 

The present paper is organized as follows. In Section 2 we review some basic 
facts about minimal formulas and give preliminary lemmas for further arguments 
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2 dim^fc (Q) - 1 if k is even and £ X, 
2 dim VI (CI) otherwise. 



(1.3) 
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in the following sections. In Section 3 by bringing the Mysovskikh theorem and a 
generalization of the Larman-Rogers-Seidel theorem together, we prove that for any 
spherically symmetric integral, the points of a minimal formula of degree 4fc + 1 on 
a particular concentric sphere have inner products all of which are rational numbers. 
In Section 4 we discuss applications of the new theory to the existence problem of 
minimal formulas: We prove that for any d > 3, there exist no d-dimensional minimal 
formulas of degree 13 and 21 for some particular integral. 

2. Preliminaries. For with < p < q < oo, let f2 = {x £ R d \ p < 

\\x\\ < q} and W be a positive weight function on f2 which is a function of ||x||. Under 
the assumption that polynomials up to a sufficient large degree are integrable with 
respect to fl for W, we consider the integral defined by 

/ f(x)W(\\x\\)dx = f 9 ( f f(rx)dp] r^W^dr, (2.1) 

where p is the surface measure on the d-dimensional unit sphere S' rf_1 . This is called a 
spherically symmetric integral. Several classical integrals like, Gaussian integral on the 
entire space, Jacobi integral on the ball, belong to the class of spherically symmetric 
integrals. Cubature formulas for spherically symmetric integrals find applications in 
various fields; for instance Lyons and Victoir jTS] recently considered cubature formula 
on Wiener space for application to mathematical finance, where cubature for Gaussian 
integral plays a key role. Hereafter we only consider spherically symmetric integrals. 

Let I be a nonnegative integer. Let Hom;(R d ) be the space of all homogeneous 
polynomials of degree I and rf — dim(Hom;(]R d )). We consider an orthonormal basis 
of Hom ; (R d ), say {P Li | 1 < % < rf}, and let P ; = (P u , . . . , P Lr <i). The bilinear form 

k 

K 2k {x,y)= Y,F 2m (x)Fj m (y). 

m=0 

is called the reproducing kernel of V^i^. 11 ) with respect to f2 for W. We note that any 
finite dimensional Hilbert space has a unique reproducing kernel; for instance see |10) . 

There is a close relationship between minimal formulas of odd degree and repro- 
ducing kernels, as the following famous theorem by Mysovskikh shows [19]; we write 
the statement only for the degree 4fc + 1 case. 

Theorem 2.1 (|19]). Assume there exists a minimal formula X of degree 4fc + 1 
for a spherically symmetric integral. Then, the formula is minimal if and only if the 
following hold: 

(i) Vx, y e X, x ±y, K 2k {x,y) = 0. 

(ii) Vx € X\{0}, w(x) = l/(2K 2k (x,x)) and w(0) = l/if 2 fc(0,0). 

We also refer the readers to [31] for relationship between minimal formulas and re- 
producing kernels. 

Theorem 12.11 includes some proven facts that the points of a minimal formula of 
degree 4fc + 1 have the antipodality and contain the origin; see [17] for details. The 
following result mentions more on the structure of minimal formulas. 

Theorem 2.2 (|12j). Assume there exists a minimal formula of degree Ak + 1 
for a spherically symmetric integral. Then the following hold: 

(i) The points are distributed over k distinct concentric spheres and the origin. 

(ii) The weights take a constant on each concentric sphere. 

(iii) The points on each concentric sphere are similar to a spherical (2fc+3) -design. 
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See also [27J for the structure of two-dimensional minimal formulas of degrees 4fc + 1 . 

To calculate the reproducing kernel of , P 2fe (R c! ), we have to compute an orthonor- 
mal basis of V$ k (M. d ) with respect to Q for W. Let Harm(R d ) be the set of all 
harmonic polynomials. Let / be a nonnegative integer. We define Harm;(M d ) = 
Harm(R d ) n Hom;(R d ). Let (fin, i = 1, . . . , dim(Harni;(R d )) be an orthonormal basis 
of Harm; (R d ) such that 

(/)i 1 j 1 (x)(fii 2i i 2 (x)dp = Si u i a 6 iuia . 



is*- 1 

It is well known (see, e.g., [21 [25]) that 

dim(Harm ! (R d )) 

£ M (x)0 M ( y ) = ^±^C«^ 2)/2) ((x, y >), V X ,yES d -i 
where C ; is the Gegenbauer polynomial of parameter A, which is explicitly given by 

rW M = (-1)' r(A + i/2)r(/ + 2A) (i - t^ /2 - x I 2 , +A _ 1/2 

1 w 2 l L(2A)L(/ + A + l/2) I! dt l[ ' 

Let k, m be nonnegative integers with k > m. We consider the polynomial space 
^2t-2mW with respect to the interval [p, q) for the weight function y~ 1 r d + im ~ 1 W(r), 
where V = f q r 4 ' 1 W(r)dr. Since T^ k _ 2m (R) consists only of even polynomials, by 
using Gram-Schmidt's orthonormalization method, we can construct an orthonormal 
basis g2m,j( r )■> j = 0, . . . ,k — m, where g 2 m,j(f ) is a polynomial of degree j in r 2 , 
such that 

6j, y = ^ [ 92 m Ar 2 )92 m ,Ar 2 )r d+im - 1 W{r)dr. 



By recalling Eq. (|2.ip . we obtain 

(Vol^-^-Va^^^^^dljf) . (y |5 d - 1 |)- 1 / 2 ^ m2 , i2 ( a ;) 52m2 , J - a (|| a ;|| 2 )^(|| a ;||)d a ; 
= |5d-i| J sd _ 1 4>2m u i 1 (x)(fi 2m2A2 (x)dp ■ — J g2 mi ,j 1 (r 2 )g2 m2 ,j 2 {r 2 )r d+im ~ 1 W(r)dr 

&ix ,12 &ji ,ji ^mi,m 2 • 

Hence (Vo\S d ~ 1 \)~ 1 ^ 2 (fi2m,i(x)g 2m ,j(\\x\\ 2 ) form an orthonormal basis of T'l k (JL d ) with 
respect to fi for W; see [51 [L2"]. 

Lemma 2.3 ([12 ]). Let d > 3. TTie reproducing kernel o/"P 2/c (IR d ) with respect to 
f2 /or W is given as follows: 



K 2 k(x,y) 



^ E E ^(ikii 2 )^(bii 2 )(ikiiiiwii) 2m d+ d 4 - 2 2 c2 - 2)/2 KhM\) ^ eRd \w> 



u| 1 m=0 j=0 

^^E^(NI 2 )WIM| 2 ) */^ = 0. 
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The following lemma is useful for calculating reproducing kernels in Section SJ 
Lemma 2.4. Let K,\ be the reproducing kernel o/'P| fc (IR.) with respect to [l,oo) 
for a weight function fi and IC2 be the reproducing kernel of V^k (K) with respect to 
(0, 1] for a weight function 7(2/) = y~ Ak ~ 2 / y) ■ Then for any x, y € (0,1] we have 



JC 2 (x,y) = (xy) 2 *^ 
or equivalently, for x, y G [1, 00), 



Kx{x,y) = (xy) 2k K 



1 1 

x 1 y 

1 1 

x' y 



Proof. Let Rj(x 2 ),j = 0, . . . , k be polynomials of degree j in x 2 , which form an 
orthonormal basis of P^fcC"^) w ith respect to [1, 00) for ji. Then by changing variable 
x — 1/y, we get 

S J , r =fR J (x 2 )R A x 2 Mx)dx = £ (y 2 % (£)) (>^, (1)) ^d y . 

Let = y 2k Rj(l/y 2 ). Then G PjfeWi an d bom the assumption 7(2/) = 

y~ 2 ~ ik [i(l/y), J Q Pj(y)Pj'(y)j(y)dy — 5jj>, i.e., Pj are an orthonormal basis of 
"P^ki^-) w ith respect to (0, 1] for 7 and hence 

' 1 \ / 1 



Since x,y G (0,1], we have l/x 2 ,l/y 2 G [l,oo). Hence from the fact Rj(l/x 2 ) 
are an orthonormal basis of V%j,(M.) with respect to [l,oo) for fi, it follows that 



Ylj=oRj(l/ x2 )Rj(l/y 2 ) is a reproducing kernel of T'JfcW with respect to [l,oo) for 
/i. The result follows by the uniqueness of reproducing kernels of Hilbert space. □ 

3. Characterizing the structure of minimal formulas. In this section we 
develop a theory of minimal formulas of degree 4/c + 1 by unifying the Mysovskikh 
theorem and a generalization of the Larman-Rogers-Seidel theorem ,22] . 

First we introduce the Larman-Rogers-Seidel theorem [M] . 

Theorem 3.1. Let X be a finite set in M. d with only two Euclidean distances 
ai,ct2 (a\ < ai) between any two distinct points. If \X\ is greater than 2d + 3, then 
a\ja\ = (m — l)/m where m is a positive integer bounded above by 1/2 + \J d/2. 

Nozaki [22] extended this theorem for points with more Euclidean distances. 

Theorem 3.2. Let X be a finite set in M. d with only s Euclidean distances 
ai, . . . , a s between any two distinct points. If \X\ > 2( <i ^] L 1 ) + 2{ d+ s s l 2 2 ^, then for 
each i = 1, . . . s, 



n 



is an integer, whose absolute value is bounded by some function of d and s. 

Now, let r%, . . . , r p be positive real numbers. Let Y be a set in M. d supported by 
Sf~ x , • • • , Sf -1 , where S~7 l is the (d— l)-dimensional sphere of radius r$. We denote 
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the ith layer of Y by Y h namely, Yi = Y H S^ 1 . Let A(YA = {(x,y)/r 2 \ x, y G 
Yi, x ^ y}. 

We can prove various theorems similar to Theorem 13.21 for antipodal finite sets 
which are supported by a sphere. 

Theorem 3.3 ([22]). Let s be an integer with s > A. Let Y be an antipodal 
finite set in Sf~ l with s inner products between any two distinct points. If \Y\ > 
4( d "^~ 3 ) + 2, then every inner product is rational. 

Theorem 13 .31 leads to the following proposition which is closely related to minimal 
formulas of degree Ak + 1. 

PROPOSITION 3.4. Let p, d be integers such that p > 2, d > Ap 2 — 2p + 1. Let Y 
be an antipodal set supported by p concentric spheres with \Y\ = 2 X)f=o ( d+ 2p-2i~ 2 
and for each i, \A(Yi)\ < 2p + 1. Then there exists Yi such that for any a € A{Y{), a 
is a rational number. 

Proof. By using the pigeonhole principle, we can choose Y\ such that 



2 ( d + 2 P - 1 - 2i 



By the assumptions of p and d, we obtain 

2 ^ (d + 2p - 1 - 2i\ 2({d + 2p-l\ fd+2p~3 



E 



, , . > , , , 

pf^\ 2p-2i j ~ p\\ 2p J \ 2p-2 



>4( d : 2P 7 2 ) + 2. (3.1) 



2p- 1 

If \A{Y{)\ < 2p, then \Y t \ < 2( d + 2p ^ 2 ) [8 . This contradicts (03]), and hence \A(Y t )\ = 
2p + 1. Thus the result follows from Theorem 13.31 □ 

The following theorem presents a number-theoretic characterization of the struc- 
ture of minimal formulas of degree 4fc + 1. 

Theorem 3.5. Let k,d be integers such that k > 2,d > Ak 2 — 2k + 1. Assume 
there exists a d-dimensional minimal formula X of degree Ak + 1 for a spherically 
symmetric integral. Then there exists a layer Xi such that every a £ A(X{) is a 
rational number. 

Proof. By Theorem 12.21 the points X are distributed over k distinct concentric 
spheres and the origin. By (|1 .3[) and Lemma 1.7], we have 

1=0 v 

By Lemma 2.3 in [T^], the reproducing kernel of V^i^) is a polynomial of degree 
k in variables ||cc|| 2 , ||y|| 2 , (x,y) 2 . If x,y G X, x ^ ±y, lie on the same layer, then 
by Theorem 12.11 K2k{x,y) is regarded as a polynomial of degree at most k in one 
variable (x,y) 2 . This implies |vl(A^)| < 2k + 1 for each layer JQ. Hence the result 
follows by Proposition [331 D 

Remark. Bannai and Damerell [S] proved that there exists no tight spherical t- 
design on S fd_1 for t > 8 and d > A except for some examples. A key idea of their 
proof is to note that the inner products of points of a tight spherical design form 
the zeros of a certain orthogonal polynomial, and moreover to show the rationality of 
these zeros by the theory of association scheme. In this sense Theorem 13.51 is similar 
to Bannai and Damerell's approach for Euclidean design. 

The usefulness of Theorem 13.51 will be clear in Section SJ 
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4. Applying the basic theory. Here we use the same notations as in the 
previous sections. We consider the following integral: 



/(•'■jn-(|k||)(b-=/ /W gE ^ 

{rreR d |l<||:r||<oo} IFII 



^ f(rx)df!) Vr 2 - 1 r-< 4fc+3 Mr. (4.1) 



This can be seen as a multidimensional generalization of the two-dimensional inte- 
gral Xu [3U] focused on, though, of course, many other generalizations can also be 
considered. The purpose of this section is to prove the following theorem. 

Theorem 4.1. There exists no d- dimensional minimal formula of degree 13 and 
21 for the integral for any d > 2. 

To prove Theorem 14.11 we use some lemmas. 

Lemma 4.2. Let d > 3. The reproducing kernel o/P^fcC^ ) with respect to \4- 1\ ) 
is given as follows: 

for any x,yeR d \{0}. 

Proof. We use the same notation g-2m,j as m Section[5] We consider the polynomial 
space f 2 *)!-2ra(^) with respect to the interval (0, 1] for the weight function y/l — r 2 . By 
the orthogonality of the Gegenbauer polynomial, 2^JV§JitC^ \r) is an orthonormal 
basis of PJ|._2 m (K) with respect to (0, 1] for \/\ — r 2 . Applying Lemma WM. we have 

ES™ <?2 m ,i(* 2 )W* 2 ) = ^- 1 K)(^) 2(fc - m) E,to m ^ 1) (l/ S )C«(lA). Hence the 
desired result follows by Lemma \2. 3 1 □ 

The following lemma determines the layers of minimal formulas for the integral 

im 

Lemma 4.3. Assume there exists a minimal formula of degree 4fc + 1 for the 
integral \4- Then the points are distributed over the origin and k concentric spheres 
of radii (cos(i7r/(2fc + 1))) _1 and the sums of weights on each concentric spheres are 
ir(2k + l)- 1 |5 d - 1 |(sin(i7r/(2fc + l))) 2 (cos(i7r/(2fc + l))) 4fe , where i = 1, • • • ,k. 

Proof. We note that for each £ = 0, • • • ,2k, 

oo 



x|| 2 ^(||a;||)da; = / dp / r 2e y/r 2 - 1 r- (4fe+3) dr 
>» Js- 1 - 1 Jl 

is^i r 1 ft + i^ 2k -"- 1 



L^) ^ At < ,4 ' 2) 



where the last equality follows by standard calculations including changing variable 
r — > y/2/(t + 1). By Theorem 12.21 the points X are distributed over the origin and 
k concentric spheres S^ -1 , • • • , S^ 1 . Let A, = J2 xeX nS d ^ 1 w ( x ) an d Ri = T 1 ■ Then 
the leftmost term of Eq. (|4.2[) can be represented as Ej=i A-i-Rf. Thus 

't+l\ 2fc_1 ~' /l -t ,, A/ 4 2k \ / „_ 1N 2fc-l-i 
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By letting m = 2k - 1 - I, Ri = (2 - R i )Rr 1 ,A i = A\S d - 1 \- 1 A i R 2k , this formula is, 
equivalently, 



2 



, -_di = ^Aj^— j , m = 0,...,2fc-l. 



This is the Gaussian quadrature formula of degree 2k— 1 for the integral tt 1 J /(i)i 

It is well known (see, e.g., |13| ) that the points and weights in the Gaussian quadrature 
formula of degree 2k — 1 are given explicitly for the univariate integral as follows: 

- 2in , ~ 4?r / . ~ 
Ri = cos — and A, = — sin ■ 



2k + 1 2k + 1 V 2fc + 1 

where i = 1, . . . , k. By using the half angle formula of the cosine function, we obtain 
Ri = r\ = (cos(«7r/(2/c + 1)))~ 2 . Moreover, the sums of weights on each concentric 
spheres are 



Ai = -R, ZK A l = — ! sin cos . □ 

4 1 2k + 1 V 2k + 1 J \ 2k + lJ 

Remark. It is basically possible to calculate the radii of layers of minimal formulas 
of degree 4fc + 1, by taking y — in Theorem 12.11 (i) and solving the equation of 
the norm ||x||. However for ordinary integrals, it seems to be rare to get "compact 
forms" of the radii of layers even for small values of k. Here "compact forms" mean 
that the radii take very complicated forms involving imaginary numbers and power 
roots, which are not desirable when we further discuss the nonexistence of minimal 
formulas: For example, it is useful to know the inner products of points in minimal 
formulas when discussing the nonexistence of minimal formulas. To do so, finding 
compact forms of radii is useful as well; see [H [T2J [57] for details. From the proof of 
Theorem 14.11 below, the readers will also see how useful such compact forms are. 

Proof of Theorem \4-l\ Assume there exists a c?-dimensional minimal formula X of 
degree ik+1 for the integral (|4.1|) . Let x,y € AT/ with x ^ ±y and = (x,y) 2 /Rj. 

By Theorem 12. 1[ Lemma 14.21 and the half angle formula of the cosine function, for 
each d > 3, we have 



0. £ i±%=lcS? (J ' + + '» ) U%) . (4.3) 

i+j=k ^ - J \ J 

i,j>0 

We regard Eq. (l4.3p as an equation of cos(2/7r/(2fc + 1)). 

(i) The Case of k = 3. We assume d > 31. According to Theorem 13.51 we can 
choose some layer of the formula, say Xi, over which the points have rational distances. 
Let f(x) = x 3 + T^x 2 — T^x — g € Q[x] be the minimal polynomial of cos(2Z7r/7) [29] . 
When dividing the right hand side of Eq. (|4.3|) by /(cos(2/7r/7)), the remainder is 

1 /2lir\ 2 

— d{d + 2){(d + 4)(d + 6)A 2 ^ m - 6(d + 4)A^ + 3} cos (— J 

+ ^d{(d + 2)(d + 4)(d + 6)A 2 ^ - 6(d + 2){d + 3)A^ + 3d} cos (^) 
+ 45(d + 2)(d 2 + 6d + 24) A^ - 15(d 2 +44)}. 
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By the minimality of /, we have 

=(d + 4)(d + 6)4^ - 6(d + 4)4^ + 3. 

=(d + 2)(d + 4)(d + 6)4^ - 6(d + 2)(d + 3)4^ + 3d. 

=(d + 2)(d + 4)(d + 6)(d + 8)(d + 10)4^ - 15(d + 2)(d + 4)(d + 6) 2 A^ 

+ 45(d + 2)(d 2 + 6d + 24)4^ - 15(d 2 + 44). (4.4) 

By multiplying the first equation in (|4.4j) by d + 2 and then subtracting the new equa- 
tion from the second equation in Eq. (|4.4j) . we see that A^-p^ = l/(d + 2). Substituting 
this value of A^j^ into the first equation in Eq. (l4.4p implies that d = —1,-6. This is 
a contradiction. We deal with the remaining cases 2 < d < 30 in Appendix [3] 

(ii) The Case of k — 5. We assume d > 91. According to Theorem 13.51 we can 
choose some layer of the formula, say A/, over which the points have rational distances. 
Let f(x) — x 5 + ix 4 — x 3 — |x 2 + jgX + ^ e Q[x] be the minimal polynomial of 
cos(2Ztt/11) [25]. When dividing the right hand side of Eq.flO]) by /(cos(2/7r/ll)), 
the remainder is 

/2l%\ /2Zvr\2 /2?7T\3 /2ZtT\ 4 

Co + Cl cos^— j+c 2 cos ^—j + C3 cos^— j + c 4 cos(^— J, 

where 



c o ==S rf {(<« + 2)(d + 4)(d + 6)(d + 8)(d + 10)(d + 12)(d + 14) (d + 16)(d + 18)^ 5 



371589121 K) ' M l " : " n " ~ """ ~~ ~ lun " r '" •"" " ; " '"" ! " ~ r " T \ 77T 

- 45(d + 2)(d + 4)(d + 6)(d + 8)(d + 10)(d + 12) (d + UfA^ 

+ 630(d + 2)(d + 4)(d + 6)(d + 8)(d + 10)(d 2 + 22d + 128)4^ 

- 3150(d+2)(d + 4)(d + 6)(d + 8)(d 2 + 16d + 84)yl 2 /S7 

+ 4725(d + 2) (d 2 + 6d + 32) (d 2 + 14d + 72) A^ - 945(d 4 + 10d 3 + 100d 2 + 440d + 4384)}, 

4 

10321920~ LV ~ ' ' ' ~ /v ~ ' ' " /v ~ ' ' ^"VHF 

- 28(d + 2)(d + 4)(d + 6)(d + 8)(d + 10)(d + 14)4^ + 210(d + 2)(d + 4)(d + 6)(d + 8)(d + 14)4 2 

- 420(d + 2)(d 3 + 24d 2 + 164d + 312)4^ + 105(d 3 + 20d 2 + 92d + 16)}, 

- 28(d + 4)(d + 6)(d + 8)(d + 10) (d + 14)4^ + 210(d + 4)(d + 6)(d 2 + 22d + 128)4^ 



ci = 1 d{(d + 2)(d + 4)(d + 6)(d + 8)(d + 10)(d + 12) (d + 14)4 



420(d + 4)(d 2 + 20d + 132)4 ^ + 105(d 2 + 18d + 152)}, 



c 4 =-r^7Ad + 2)(d + 4){(d + 6)(d + 8)(d + 10)4 3 - 15(d + 6)(d + 8)4 2 + 45(d + 6)4^ - 15}. 



c 3 =4^80^ + 2 ){( d + 4 )( d + 6 )( d + 8 K d + 10 )^ " 15 ^ d + + 6 ^ d + ^ A VRl 
+ 45(d + 4) 2 4 v777 - 15(d+ 1)}, 
1 

46080 ' 

By the minimality of /, we have Cq = 0, c\ = 0, c 2 = 0, c 3 = 0, c 4 = 0, and therefore 
/ 3 ,4 := c 3 - c 4 =^{(^ 3 + 12rf 2 + 44d + 48)4^ - 6(d 2 + 6d + 8)4^ + 3(d + 2)}, 
h,x ~ c 2 - ci =^{(^ 3 + 12rf2 + 44d + 48)4^ - 3(2d 2 + 13d + 18) + 3(d + 3)}. 
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Since 

n . , _ d{(d + 2)A vm -l} 
- Ha ~ h,x ^ , 

we obtain A^j^ = l/(d + 2). By substituting A^rp^ into /3,4, we obtain 

_ d(d+l)(d + 6) 
/3 ' 4 "" 1536(d+2) <0 - 

This is a contradiction. We deal with the remaining cases 2 < d < 90 in Appendix [3] 
Remark. When 2fc+l is a prime number, the minimal polynomial of cos(27ri/(2fc+ 
1)) is uniquely determined, not depending on I [59]. Using this fact, we proved there 
exist no minimal formulas of degrees 13 and 21 for the integral (|4.1j) . Similar but 
more extensive considerations will enable us to discuss the nonexistence of minimal 
formulas for other degrees, though we do not pursue it here. 

5. Conclusion. Unifying the Mysovskikh theorem and a generalization of the 
Larman-Rogers-Seidel theorem, we obtained a new necessary condition for the exis- 
tence of minimal formulas of degree 4fc + 1 for spherically symmetric integrals. And 
thereby we solved the existence problem of minimal formulas of degrees 13 and 21 for 
the integral ()4.1j) . In this paper we only focused on the particular integral. We hope 
that our approach could also be applied for other spherically symmetric integrals. 

Appendix A. Small dimensional cases of the proof of Theorem 14.11 

In this section we finish off the remaining cases of the proof of Theorem 14. II The 
following lemma plays a key role. 

Lemma A.l. Let d > 3. Assume there exists a d-dimensional minimal formula 
of degree 4fc + 1 for the integral \4- Then, 

8 . (2j + l)7r\ 2 f/d+2m-l\ (d + 2m - 3\ 1 , Ai , 

is an integer. 

Proof. We use the same notations ri, S^ -1 , Ai as in Section [4] According to 
Theorem 12.21 (ii), let w\ = w(x) for every iglfl Sf _1 . By combining Theorem 12. II 
(ii) and Lemma 14.21 we have W\ = \/{2K 2 k{x, x)). Thus by Lemma[ 

2 k k—rn , , . _ , , \ \ 2 



* n *-i - -5tt(*stt) S t *^{<* , (« 5 tt)} ^-<u 

v / m— j—0 

8 y y ( ,■ (2 J + Ifr \ 2 d + 4m - 2 ((d -2)/2) n , 
~2fcTT ^ ^ ^ bln 2fc + l J d-2 C2m (1) 



(2j + l)ir\ ( fd + 2m- l\ /d+2m-3 
sin v 1)1 1 ' 



m— j—0 

q k k—m 

2fc + l^^V 2fc + 1 J IV 2m / V 2m - 2 

m— j—0 

where the last equality follows by the fact d+ ^ 2 C { ^ 2)/2) {I) = dim(Harm 2 „ l (K d )) 
(see 0). □ 
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Let n be an positive integer. For j = 0, . . . , k, we obtain approximate sine func- 
tions 




where tt± are positive real numbers with 7r_ < n < n + . These functions satisfy 
Si(n;k,d) < sin((2j + l)7r/(2fc + 1)) < S 2 (n;k,d). By using Lemma ED and the 
above equations, we have 

m=0 j=0 k v 7 v 7 

8 , f (d + 2m- l\ /d+2m-3\l 

m=0 j=0 k v 7 v 7 J 

By choosing the suitable values n and ir± , we see that (|A.1|) are not integers when (i) 
fc = 3, 3 < d < 30, (ii) k = 5, 3 < d < 90. 

For example, for (fc, d) = (3,20), by substituting n = 5, n + = 3.14160 and 
7r_ = 3.14159, we obtain 

47868.2 < N\ < \X\\ < N 2 < 47868.8. 

When (k,d) = (5,80), by substituting n = 11, tt + = 3.1415926535898 and tt_ = 
3.1415926535897, we obtain 

318122993450.96 < Ni < \Xi\ < N 2 < 3.18122993450.99. 

It remains to consider the case when d = 2. Since the radii and weights of minimal 
formulas of degree 4fc + 1 are determined by Lemma 14.31 By applying these values 
to the equation in [S] Theorem 3.1.1 (2)], we easily see that there exists no two- 
dimensional minimal formula, which completes the proof. 
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